Quantum circuit for three-qubit random states 
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We explicitly construct a quantum circuit which exactly generates random three-qubit states. 
The optimal circuit consists of three CNOT gates and fifteen single qubit elementary rotations, 
parametrized by fourteen independent angles. The explicit distribution of these angles is derived, 
showing that the joint distribution is a product of independent distributions of individual angles 
apart from four angles. 
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I. INTRODUCTION 

Quantum information science (see e.g. [l[ and refer- 
ences therein) has received an increased attention in re- 
cent years due to the understanding that it enables to 
perform procedures not possible by purely classical re- 
sources. Experimental techniques to manipulate increas- 
ingly complex quantum systems are also rapidly progress- 
ing. One of the central issues is on the one hand to con- 
trol and manipulate delicate complex quantum states in 
an efficient manner, but on the other hand at the same 
time to prevent all uncontrollable influences from the en- 
vironment. In order to tackle such problems, one has 
to understand the structure and properties of quantum 
states. This can be done either through studies of par- 
ticular states in a particular setting, or through focusing 
on the properties of the most generic states. 

Random quantum states, that is states distributed ac- 
cording to the unitarily invariant Fubini-Study measure, 
are good candidates for describing generic states. Indeed, 
they are typical in the sense that statistical properties of 
states from a given Hilbert space are well described by 
those of random quantum states. Also, they describe 
eigenstates of sufficiently complex quantum systems @ 
as well as time evolved states after sufficiently long evolu- 
tion. Not least, because random quantum states possess 
a large amount of entanglement they arc useful in cer- 
tain quantum information processes like quantum dense 
coding and remote state preparation [H, 0] . They can be 
used to produce random unitaries needed in noise esti- 
mation Q and twirling operations [||. In addition, as 
random states are closely connected to the unitarily in- 
variant Haar measure of unitary matrices, the unitary in- 
variancc makes theoretical treatment of such states sim- 
pler. 

Producing such states therefore enables to make avail- 
able a useful quantum resource, and in addition to span 
the space of quantum states in a well-defined sense. 
Therefore several works have recently explored different 
procedures to achieve this goal. It is known that generat- 
ing random states distributed according to the exact in- 
variant measure requires a number of gates exponential in 



the number of qubits. A more efficient but approximate 
way to generate random states uses pseudo-random quan- 
tum circuits in which gates are randomly drawn from a 
universal set of gates. As the number of applied gates 
increases the resulting measure gets increasingly close 
to the asymptotic invariant measure Q. Some bipartite 
properties of random states can be reproduced in a num- 
ber of steps that is smaller than exponential in the num- 
ber of qubits. Polynomial convergence bounds have been 
derived analytically for bipartite entanglement [1, d, |T(| 
for a number of pseudo-random protocols. On the nu- 
merical side, different properties of circuits generating 
random states have been studied [ill ]. In order to quan- 
tify how well a given pseudo-random scheme reproduces 
the unitarily invariant distribution, one can study aver- 
ages of low-order polynomials in matrix elements [12j | . In 
particular, one can define a state or a unitary fc-design, 
for which moments up to order k agree with the Haar dis- 
tribution [l3|, EH ■ Although exact state fc-designs can be 
built for all k (see references in [l4|) they are in general 
inefficient. In contrast, efficient approximate fc-designs 
can be constructed for arbitrary k (for the specific case 
of 2-design see @). 

The pseudo-random circuit approach can yield only 
pseudo-random states, which do not reproduce exactly 
the unitarily invariant distribution. The method has 
been shown to be useful for large number of qubits, where 
exact methods arc clearly inefficient. However, for sys- 
tems with few qubits, the question of asymptotic com- 
plexity is not relevant. It is thus of interest to study 
specifically these systems and to find the most efficient 
way - in terms of number of gates - to generate ran- 
dom states distributed according to the unitarily invari- 
ant measure. This question is not just of academic in- 
terest since, as mentioned, few-qubit random unitaries 
are needed for e.g. noise estimation or twirling opera- 
tions. Optimal circuits for small number of qubits could 
also be used as a basic building block of pseudo-random 
circuits for larger number qubits, which might lead to 
faster convergence. In addition, systems of few qubits 
are becoming available experimentally, and it is impor- 
tant to propose algorithms that could be implemented on 
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such small quantum processors, and which use as little 
quantum gates as possible. Indeed, quantum gates, and 
especially two-qubit gates, are a scarce resource in real 
systems which should be carefully optimized. 

In this paper we therefore follow a different strat- 
egy from the more generally adopted approach of us- 
ing pseudo-random circuits to generate pseudo-random 
states, and try and construct exact algorithms generat- 
ing random states for systems of three qubits. In the lan- 
guage of fc-designs such algorithms are exact oo-designs. 
We present a circuit composed of one-qubit and two- 
qubit gates which produces exact random states in an 
optimal way, in the sense of using the smallest possi- 
ble number of CNOT gates. The circuit uses in total 
3 CNOT gates and 15 one-qubit elementary rotations. 
Our procedure uses results recently obtained [l5[ which 
described optimal procedures to transform a threc-qubit 
state into another. Our circuit needs 14 random numbers 
which should be classically drawn and used as parame- 
ters for performing the one-qubit gates. The probability 
distribution of these parameters is derived, showing that 
it factorizes into a product of 10 independent distribu- 
tions of one parameter and a joint distribution of the 4 
remaining ones, each of these distributions being explic- 
itly given. Since we had to devise specific methods to 
compute these distributions, we explain the derivation 
in some details, as these methods can be useful in other 
contexts. 

After presenting the main idea of the calculation in 
Section [H] we start by treating the simple case of two- 
qubit states in Section IIIII We then turn to the three- 
qubit case and first show factorization of the probability 
distribution for a certain subset of the parameters (Sec- 
tion |IV|, the remaining parameters being treated in Sec- 
tion |Vl The full probability distribution for three qubits 
is summarized in Section IVII 



II. THE QUANTUM CIRCUIT 

Formally, a quantum state can be considered as an 
element of the complex projective space (DP^" 1 , with 
N = 2™ the Hilbcrt space dimension for n qubits [l6j |. 
The natural Riemannian metric on (DP^ -1 is the Fubini- 
Study metric, induced by the unitarily invariant Haar 
measure on U(N). It is the only metric invariant un- 
der unitary transformations. To parametrize (DP^ -1 
one needs 2N — 1 independent real parameters. Such 
paramctrizations are well-known, for instance using Hur- 
witz paramctrization of U(N) [l7[. However, they do 
not easily translate into one and two-qubit operations, 
as desired in quantum information. In Rcf. [15| . optimal 
quantum circuits transforming the three-qubit state |000) 
into an arbitrary quantum state were discussed. In the 
case of three qubits, a generic state can be parametrized 
up to a global phase by 14 parameters. The quantum 
circuit requiring the smallest amount of CNOT gates has 
three CNOTs and 15 one-qubit gates depending on 14 



independent rotation angles. From [15[ it is possible (see 
Appendix) to extract the circuit depicted in Fig. [TJ ex- 
pressed as a series of CNOT gates and single qubit ro- 
tations, where Z-rotation is Zg = exp (— ia z 9) and Y- 
rotation is Yg = exp (— icr y 0) with a yjZ the Pauli matri- 
ces. The circuit allows to go from 1 000) to any quantum 
state (up to an irrelevant global phase). It therefore pro- 
vides a paramctrization of a quantum state \ip) by angles 

01, ... , #14. 

In order to generate random vectors distributed ac- 
cording to the Fubini-Study measure, it would of course 
be possible to use e. g. Hurwitz parametrization to gen- 
erate classically a random state, and then use the proce- 
dure described in [l5[ to find out the consecutive steps 
that allow to construct this particular vector from 1 000) . 
However this procedure requires application of a specific 
algorithm for each realization of the random vector. In- 
stead, our aim here is to directly find the distribution 
of the di such that the resulting \ip) is distributed ac- 
cording to the Fubini-Study measure. This is equiva- 
lent to calculating the invariant measure associated with 
the paramctrization provided by Fig. [T] in terms of the 
angles 6\, . . . , 0m. Geometrically, the Fubini-Study dis- 
tance Dps is the angle between two normalized states, 
cos (Dps) = KV'I'WI- The metric induced by this distance 
is obtained by taking \cf)) = + \dip), getting 

2 _ (^)(d^,d^)-(^,dy>){d?M) 

fa*)* [1) 

where (, ) is the usual Hcrmitian scalar product on 
fl8| . If a state is parametrized by some parameters 
0i , 62 , . . • then the Riemannian metric tensor gij is such 
that ds 2 = gijd^id^j and the volume form at each point 
of the coordinate patch, directly giving the invariant mea- 
sure, is then given by dv = ^/det(g) FJ d0j . Thus the joint 
distribution P(6) of the 0i is simply obtained by calcu- 
lating the determinant of the metric tensor given by (fTJ) 
with the parametrization \ip) = \ip(d)}, 9 = (61, ... , #14). 
Unfortunately the calculation of such a 14 x 14 determi- 
nant for n = 3 qubits is intractable and one has to resort 
to other means. Let us first consider the easier case of 
n = 2 qubits, where by contrast the calculation can be 
performed directly. 

III. A SIMPLE EXAMPLE: THE TWO-QUBIT 
CASE 

A normalized random 2-qubit state depends, up to 
a global phase, on 6 independent real parameters. A cir- 
cuit producing \ip) from an initial state |00) is depicted 
in Fig. [2j One can easily calculate the parametrization 
of the final state \i/j(9i, . ■ ■ , 0e)) in terms of all six an- 
gles, thus directly obtaining the metric tensor from 
(fl"j). Square root of the determinant of g then gives an 
unnormalizcd probability distribution of the angles as 

P(0) = I cos 2 20i sin 20i sin 20 3 sin 20 5 1 (2) 
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FIG. 1: Circuit C for three-qubit random state generation. 
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FIG. 2: Circuit for two-qubit random state generation. 



(see also [lj|). Several observations can be made about 
this distribution. First, the three rotations applied on 
the first qubit (top wire in Fig. |2J) after the CNOT gate 
represent a random SU(2) rotation, for which the Y- 
rotation angle is distributed as P{9 3 ) ~ | sin 2^3 1 and 
the Z-rotation angles are uniformly distributed (l7j . Sec- 
ondly, angle 0i gives the eigenvalue of the reduced den- 
sity matrix, A = sin 2 9i , for which the distribution is well 
known, see, e.g., [2G|. The third observation is that, 
remarkably, the joint distribution @ of all 6 angles fac- 
torizes into 6 independent one-angle distributions. 



IV. FACTORIZATION OF THE THREE-QUBIT 
DISTRIBUTION FOR ANGLES 9 7 TO 14 



Let us denote by C the circuit of Fig. Q] and by 
C(9) the unitary operator corresponding to it, so that 
\ip(0)) = C(0)|OOO). Because circuits C span the whole 
space of 3-qubit states, any unitary 3-qubit transforma- 
tion V maps parameters 9 to new parameters 9 such 
that V\ip(9)) = \ip(9)). We denote by C the circuit 
parametrized by angles 9 corresponding to performing 
C followed by V . It is associated with the unitary oper- 
ator C{9) such that C(0)|OOO) = VC(9)\000). Unitary 
invariance of the measure implies for P(9) that 



P(9) = P(9)\J\, 



(3) 



with J the Jacobian of the transformation 9 i— > 9 and | . | 
denotes the determinant. Note that Eq. © is not a sim- 
ple change of variables, as the same function P appears 
on both sides of the equation. The Jacobian matrix J for 
transformation V from angles 9 to 9, V\ip(0)) = \i/j(9)}, 
tells how much do the angles 9 of \ip(9)) change if we 
vary angles 9 in \ip{9)) keeping transformation matrix V 
fixed. Choosing V that sets some angles 6j in circuit C 
to a fixed value, say zero, and at the same time showing 
that \J~\ depends only on these angles 6j, would prove 
factorization of P(9) with respect to angles 9j through 
Eq.H 



Let us now turn to our main issue, which is the distri- 
bution of angles in the threc-qubit case. In order to have 
an indication whether the distribution of an angle (9; fac- 
torizes, we numerically computed the determinant det(g) 
of the metric tensor as a function of 9i with the other an- 
gles fixed. We also numerically computed the marginal 
distribution of 9i by using the procedure given in the ap- 
pendix to find the angles corresponding to a sample of 
uniformly distributed random vectors. If the distribution 
for a given angle 9i factorizes, these two numerically com- 
puted functions should match (up to a constant factor). 
This is what we observed for all angles but four of them 
(angles 9 3 to 9 6 ). 

In order to turn this numerical observation into a rig- 
orous proof, we are going to show in this section that the 
distributions for angles 9j to #14 indeed factorize. In the 
next section we will complete the proof by dealing with 
the cases 9\ to 9§. The explicit analytical expression of 
the probability distribution for individual angles will be 
given in Section IVTl 



A. Gates 7-12 and 14 

The simplest case is that of gates at the end of the 
circuit C of Fig. [TJ e.g., gate Yg 14 . For V we take Y- 
rotation by angle — u on the third qubit, V = Y- u - It 
defines a mapping 9^9 such that 0j = 9i for i < 13 
and 9i4 = #14 — u. Matrix elements of the Jacobian, i.e. 
partial derivatives Jjk = 09j/d9k, are equal to the Kro- 
neckcr symbol Sjk- The Jacobian is equal to an identity 
matrix and its determinant is one. Equation taken at 
u = 0i4 then gives P(9 1 , 9 13 , 0) = P(9 U 9 13 ,6 14 ), 
from which one concludes that the distribution for 9u 
factorizes and is in fact uniform (unless noted otherwise 
P's are not normalized). The same argument holds for 
the two other rotations by angles 0i2 and 0g applied at 
the end of each qubit wire. 

Proceeding to angle 6s one could use V = Y_ Us Z_ Ug 
applied on the first qubit and show that the Jacobian de- 
pends only on 0g and 9g, while at u% = 9$ and ug = 9g one 
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gets #8 = Og = 0, from which factorization of 6*§ would fol- 
low from Eq. [3J There is however a simpler way. Observe 
that the three single-qubit gates with angles 07, 6% and 
0g on the first qubit span the whole SU(2) group. There- 
fore, for any one-qubit unitary V, gates VZg g Yg s Zg 7 can 
be rewritten as Zg Y@ Z§ , without affecting other 0's. 
The distribution of these three angles must therefore be 
the same as the distribution of corresponding SU(2) pa- 
rameters. Note that the same argument can be applied 
in the 2-qubit case of Section Hm As a consequence, the 
distribution of angles for gates Z — Y — Z at the end 
of the circuit should be the same in both cases, that is 
the distribution of 6j is uniform while that of 6s is pro- 
portional to | sin 26*8 1 . Similarly, one can show that the 
distribution for the angles 0io to #12 is the same as for 
angles 07 to 0g. 



B. Gate 13 

As opposed to gates 10-12, for gate 13 we can not 
use the analogy with the 2-qubit circuit (Figf5]) because 
the two gates 13 and 14 on the third qubit do not span 
the whole SU(2) group. Therefore a different argument 
should be used. In what follows we show that the joint 
distribution for angles #13 and 0i 4 can be factorized out 
of the full distribution. Since it has been shown in Sub- 
section |TV3] that angle 6*14 factorizes, this will prove that 
the distribution for #13 also factorizes. 

Using V = Z— Ul3 Y— Uli on the third qubit we can set 
#13 and #14 to zero with the choice U13 = #13 and Uu = 
6>i 4 . Our goal is to show that \J\ depends only on #13 
and 6*14. We can formally consider each angle 6i as being 
a function 0»(0; U13, Uu) of the initial 9 as well as of 
the parameters Ui 3 ,i 4 through C(0)|OOO) = UC*(0)|OOO). 
To calculate matrix elements of J for our choice of V 
evaluated at 1*13 = #13 and U14 = #14, we must obtain 
the first-order expansion in e of the quantities 



, 0fc-i, Qk + e, 0fc+i, 



'14, C13, "14 



)■ (4) 



Some angles 0j are very simple. We immediately see that 
when varying angles 0fc<i 2 , that is taking k < 12 in Eq[41 
angles 0j<i2 do not change (i.e 6j = 6j). The correspond- 
ing 12 x 12-dimcnsional subblock in J is therefore equal 
to an identity matrix. Similarly, taking k — 13 in Eqf4] 
we see that 0j<i3 do not change. The corresponding col- 
umn in J is therefore zero apart from 1 on the diagonal. 
The Jacobian thus has a block structure of the form 



\J\ 



1 B 
A 



\A\, 



(5) 



where 1 is a 13 x 13-dimensional identity matrix and 
A is a 1 x 1-dimensional block with partial derivative 
8614/8614. The angle 614 given by Eq.(Tj||) is obtained by 
varying angle #14 by e. The condition that C(#)|000) = 
VC{8)\000) iBY SiA Z Sia \x) = Z- 9l3 Y e Zg 13 \x), where \x) = 



CNOT 23 |^} is a state after the third CNOT acts on 
(counting from the left in Fig. [T]), explicitly given by 



cos 0i 1 00a) + sin 0i 1 1) (sin20 2 |O) + cos20 2 |l)) 



with 



|a) 



„2i0i 



e 2W(> (cos 3 sin 5 + e 2i 



cos ft 3 cos f 5 — e""" sin ft 3 sin 5 ) |0) 

sin 3 cos 05) |1) 
(sin 3 cos 5 - e 2Wl cos 3 sin 5 ) |0) 
s 2ie 6 t 2ie 4 cog da CQS g 5 + gin <9 3 sin 6> 5 ) ] 1> . 



(6) 
(7) 
(8) 



\x) is therefore determined by a, (3 and (3, where \[3) — 

a-x\(3). Similarly, \x) is determined by a,/3 and (3. Pro- 
jecting these conditions for \(3) and \f3) on the computa- 
tional basis {|0), |1)}, eliminating unwanted variables, we 
get the following two equations: 



cos e — k sin e e 
cos e + k sin e e _ 



cos0i 4 e l6 ' 13 + k sin 0i4e — 1 13 



(9) 

one with k = 1, one with k = — 1. Expanding these equa- 
tions to first order in e yields 9 14 = e cos 20i3 + o(e 2 ) and 
013 = + o(e 2 ). The derivative 8614/8614 is therefore 
equal to cos 20i3, which completes the proof. Inciden- 
tally, we also see that the distribution of 0i3 is propor- 
tional to I cos 2013 1 . 



V. JOINT THREE-QUBIT PROBABILITY 
DISTRIBUTION FOR ANGLES 0% TO 6 

In the preceding section, we have shown that the distri- 
bution for angles 6j to 0i4 factorizes. As was mentioned, 
numerical observations indicated us that the distribution 
for angles 0i and 2 should also factorize, but that it is 
not the case for the joint distribution of 03, . . . , 6$. 

As we were not able to directly prove by the same 
methods as above that the distributions for 0i and 2 
factorize, we use a different strategy. Namely, we first 
assume that this factorization is true, then we compute 
the distributions under this assumption, and the knowl- 
edge of the answer allows us to prove a posteriori that it 
is indeed the correct probability distribution. 

If the factorization holds, the distribution for 0i and 
2 is easily calculated from the matrix g using symbolic 
manipulation software, by replacing angles 6j, j > 3, in 
g by suitably chosen simple values, so that the 14 x 14 
determinant giving the volume form can now be handled. 
This yields, up to a normalization constant, 



Pi(0i) = cos 5 fti sin" pi 



5 p „:„9 , 
5 o/i„ c,;v,3 , 



P 2 (0 2 ) = cos b 20 2 sin J 20 2 . 



(10) 
(11) 



The joint distribution of 03, . . . , 0g can not be further 
factorized, and requires heavy calculations. Indeed, even 
replacing all angles but 03, 0g by numerical values the 
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determinant det(g) of the metric tensor given by {T]) still 
depends on 4 variables, which is too much for it to be 
evaluated by standard software. We thus proceed as fol- 
lows. First one can show that det(<?) can be put under 
the form 

10 6 8 6 

det(g) = a pqrs cos(2p0 3 + 2q9 4 

p— — 10 q— — 6 r— — 8 s— — 6 

+2r0 5 + 2s0 6 ), (12) 

with the sums running over all q, r but only even val- 
ues of p and s. Because of the parity of cos, there are 
M = 8509 independent coefficients a pqrs . Evaluating nu- 
merically the determinant at M random values of the an- 
gles one gets anMxM linear system that can be solved 
numerically. If the values of the coefficients of the matrix 
gij are multiplied by a factor 4, then one is ensured (from 
inspection of det(g)) that the a pqrs are rationals of the 
form fc/2 9 , k £ Z. This allows to deduce their exact value 
from the numerical result. We are left with 6998 nonzero 
terms in det(g), and terms with odd q or r do not exist. 
We then suppose that yMetQ?) can be expanded as 

v/det(^)= (13) 

5 3 4 3 

y^ y^ b pqrs e^ 2pe3+2qf>4+2r95+2se6 ' ) 

p— — 5 q— — 3 r— — 4 s=— 3 

This assumption is validated a posteriori, since a solution 
of the form (fl~3|) can indeed be found. There are 4851 co- 
efficients b pqrs , which can be obtained by identifying term 
by term coefficients in the expansion of (y/det(g)) 2 and 
det(g). We have to solve a system of quadratic equations 

010,6,8,6 = 65343 (14) 

<2l0,6,8,5 = 6534265343 + 6534365342 

010,6,8,4 = 6534165343 + ^5342 ^5342 + 6534365341 

The first equation is quadratic and fixes an overall 
sign. Equation k + 1 is linear once the values ob- 
tained from equations 1 to k are plugged into it. Start- 
ing with the highest-degree term (p,q,r,s) = (5,3,4,3) 
one can thus recursively solve all equations. There 
are only 1320 non-zero coefficients b pqrs . Gather- 
ing together terms (±p, ±g, ±r, is) one can simplify 
the sum (IT3l) to a sum of 96 terms of the form 
Cp q rs cos{p03) cos(g#4) cos(r0s) sin(s^e)' Expanding this 
expression in powers of cos(205) and sin(205) and simpli- 
fying separately each coefficient we finally get 

P(0 3 , 4 , 5 , 6 ) = sin 20 5 sin 40 3 sin 2 <pi cos ip 2 , (15) 

where (a\/3) = cos^i and ((3\f3) = cos ip 2 with \a), \(3) 
given by Eq. ([7]). Recall that |/3) is the bit-flip transform 
of |/3), |/3) = cr x Note that sin20 3 = (a\{3). Angles 
Lpifl can be obtained from 

COS 2 (fii = (C 4 C 5 C 6 - S 4 S 6 ) 2 + (C 3 C 4 S 6 + C 3 S4C 5 C 6 ) 2 
cosip 2 = -s 3 s 4 s 6 + C 6 (s 3 C 4 C5 - C3S5), (16) 



where c, = cos20i and s, = sin20i. We do not have a 
general argument to explain this remarkable expression 
of the distribution in terms of the scalar products of |a), 
|/3) and |/3). 

To complete the proof for the joint distribution 
P(0i,...,06) it remains to be checked that the de- 
terminant of the metric tensor g with angles 07 to 
014 replaced by constants is indeed proportional to 
Pi(0i)P2(02)P(03, 04,05,0e)- This a posteriori verifica- 
tion is easier to handle symbolically than the full a pri- 
ori calculation of the 14 x 14 determinant. Indeed, the 
determinant can first be reduced to an 8 x 8 determi- 
nant by Gauss- Jordan elimination. The remaining de- 
terminant can be expanded as a trigonometric polyno- 
mial. Although symbolic manipulation software do not 
allow to simplify the coefficients of this polynomial, they 
are able to check that these coefficients match those of 
the expected distribution. We proved in that way that 
the difference between the determinant det(gr) and our 
expression is identically zero. This gives a computer- 
assisted but rigorous proof for the distribution of angles 

01 to 6 . 

VI. TOTAL THREE-QUBIT PROBABILITY 
DISTRIBUTION FUNCTION 

Gathering together the results of the previous sections 
we obtain that the joint distribution P(8) can be factor- 
ized as 

14 

P(0) = |Pl(0l)P 2 (0 2 )P(03, 04, 05, 06) 11^(^)1- ( 1? ) 

i=7 

The joint distribution P(03, 04, 05, 06) has been derived 
in the previous section and is given by Eq. (|15[) . The dis- 
tribution for 0i and 9 2 is given by Eqs. pJ3l) - (fTT|l . Given 
the factorization l|17p. it is easy to calculate the remain- 
ing Pj(0i) for each i = 7, . . ., 14 as was done for 0i and 

02 in the previous section: replacing angles 0j, j ^ i, in 
g by suitably chosen simple values, the 14 x 14 deter- 
minant giving the volume form can be easily evaluated 
by standard symbolic manipulation. This yields, up to a 
normalization constant, 

14 

~[p i (9 l ) = sin 20 8 sin 20ii cos 20i3. (18) 

»=7 

The knowledge of the angle distribution (jTTJ) allows to 
easily generate random three-qubit vectors using the cir- 
cuit of Fig. [TJ Angles 0i, 9 2 and 07 to 0i4 can be drawn 
classically according to their individual probability dis- 
tribution. Angles 03 . . . , 06 can be obtained classically 
from the joint distribution (jT5|) by, for instance, Monte- 
Carlo rejection method (that is, drawing angles 03 to 06 
and a parameter x £ [0,p] at random, and keeping them 
if P(0 3 , 04,05, 6 ) < x). Bounding P(0 3 , 4 , 5 , 6 ) from 
above by p = 0.85 yields a success rate of about 12%. 
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VII. CONCLUSION 

In this work, we constructed a quantum circuit for gen- 
crating three-qubit states distributed according to the 
unitarily invariant measure. The construction is exact 
and optimal in the sense of having the smallest possible 
number of CNOT gates. The procedure requires a set 
of 14 random numbers classically drawn, which will be 
the angles of the one-qubit rotations, and whose distri- 
bution has been explicitly given. Remarkably, we have 
shown that the distribution of angles factorizes, apart 
from that of four angles. The circuit can be used as a 
three-qubit random state generator, thus producing at 
will typical states on three qubits. It could be also used 
as a building block for pseudo-random circuits in order to 
produce pseudo-random quantum states on an arbitrary 
number of qubits. At last, it gives an example of a quan- 
tum algorithm producing interesting results which could 
be implemented on a few-qubit platform, using only 18 
quantum gates, of which 15 are one-qubit elementary ro- 
tations much less demanding experimentally. 



where are one-qubit states, \uii) is a one-qubit state 
orthogonal to \u>\) and |£)23 is a two-qubit state of the 
second and third qubits. The angle 9g is chosen such 
that the Z-rotation of angle —9g eliminates a relative 
phase between the coefficients of the expansion of 
into |0) and |1). (Note that because we are using the 
circuit in the reverse direction the angles of rotations 
have opposite signs). A subsequent Y-rotation with an- 
gle —6*8 results in the transformation \oji) — ► |0) (up to 
a global phase). Similarly, rotations of angles —#12 and 
—9n rotate \u>2) into |0). After applying rotations of an- 
gles —6*8, —6*g, —6*ii and — 0\2 the state has become of 
the form |V>') = IOO7) + |l)(|07i) + 1-72)) (up to nor- 
malization). Two rotations on the third qubit of angles 
—#13 and —0i4 are now chosen so as to rotate I71} into 
some new state I7'} while I72) is rotated, up to nor- 
malization, into |t'') . It was shown in [151 ] that this 
can always be done by writing the normalized (71,2) as 
|7i,2) = cos^i^lO) + e 1 ^ 1 ' 2 sin0i j2 |l), and then #14 is a 
solution of 



- tan (26>u) = - 



cos 20i + cos 2</>2 



sin 20i cos £1 + sin 20 2 cos £2 



(A.2) 
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APPENDIX: THE PARAMETRIZATION 
CORRESPONDING TO THE THREE-QUBIT 
CIRCUIT 

In this Appendix, we explain how to obtain the angles 
6i of the circuit (Fig. [1]) for a given based on the 
discussion in [l^]. This justifies the use of these angles as 
a parametrization of the quantum states. We start from a 
state \ip), and transform it by the inverse of the different 
gates of Fig. [1] to end up with 1 000) , specifying how the 
angles 9i are obtained in turn. More details can be found 
in [15| . A generic three-qubit state can be written in a 
canonical form as a sum of two (not normalized) product 
terms (2l| . 



\il>) = \ujiuj 2 uj 3 ) + \uJi)\£,) 



2:S- 



(A.l) 



while 6*i3 = — (61 + 62)/ '4, where <5's are relative phases 
in ^-e 14 |7i,2) = c 1 * 1 ' 2 cos k|0) +sinK|l). Acting with a 
CNOT23 gate on the resulting state one obtains a quan- 
tum state for the three qubits of the form \ip") = |00xi) + 
IIW4X2), with \2 = 7'- The Z-rotation angle —Q\q on the 
second qubit is now determined so as to eliminate a rel- 
ative phase between the expansion coefficients of IW4), 
making them real up to a global phase. On the third 
qubit we now apply three rotations of angles —9 4, —#5, 
and — 9q to bring \xi) to \x') and ^2) into <t x |x')j elimi- 
nating also a relative phase. Then a CNOT13 is applied. 
At this point (after the second CNOT in Fig. [TJ counting 
from right, but without the #7 rotation), the state has be- 
come of the form \tp"') = cos 9\ lOOa^} + e lr sin 9\ \\u)§lo§), 
where the one-qubit states \lu^} and \loq) are normalized 
and real. With 9j wc now eliminate the relative phase 
t, and with an Y-rotation of angle —63 the third qubit 
is brought to the state |0). Then the combination of two 
Y-rotation of angles 6* 2 and —82 with a CNOT12 brings 
the second qubit to |0), and the last rotation of angle —9\ 
on the first qubit yields the final state |000) . Note that in 
the circuit of Fig. Q] the two Z-rotations of angles #7 and 
6>io commute with CNOT gates if they act on the control 
qubit. This is the reason why the rotation of angle 9y can 
be applied at any point between 9\ and 9$ and, similarly, 
#10 can be applied at any point between #2 and 9\\. 
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